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We present an analytical solution for the tunneling process in a piecewise linear and quadratic 
potential which does not make use of the thin-wall approximation. A quadratic potential allows 
for smooth attachment of various slopes exiting into the final minimum of a realistic potential. 
Our tunneling solution thus serves as a realistic approximation to situations such as populating a 
landscape of slow-roll inflationary regions by tunneling, and it is valid for all regimes of the barrier 
parameters. We shortly comment on the inclusion of gravity. 



Vacuum decay is one of the most drastic environmental 
changes in field theory, which often proceeds via nucleat- 
ing a bubble of true vacuum inside a sea of false vacuum 
by a quantum mechanical tunneling event. Tunneling 
processes (in first order phase transitions) play a vital 
role in many aspects of high-energy theory and cosmol- 
ogy. The groundwork for the computation of tunneling 
amplitudes was laid many years ago by Coleman and de 
Luccia (CdL) [TJ [2]. Their method uses the Euclidean 
path integral to calculate what is now known as the CdL 
instanton, often appearing to be the stationary point of 
minimal action. The CdL analysis consist of a single real 
scalar field in a potential with a false and a true vacuum 
located at the position of the corresponding minima of 
the scalar potential, <fi + and 0_ respectively, see Figure]!] 
The tunneling probability per unit volume T/V = Ae~ , 
can be conveniently expressed |12j in terms of the Eu- 
clidean action Se of the 0(4) symmetric so-called bounce 
solution (/"Bounce as B = S E [(j) Bounce } - S E [c/) + } . In the 
interior of the nucleating bubble, the scalar field exits 
at some point on the slope towards the minimum, while 
outside of the bubble, beyond the bubble wall, it still sits 
in the false minimum. 

In general, the equation of motion for the bounce is dif- 
ficult to solve analytically. This lead to the development 
of the thin- wall approximation pQ. In this limit of small 
energy difference between the true and false vacuum, the 
bubble wall becomes infinitely thin, and an approximate 
solution can be found. 

To the best of our knowledge, an exact tunneling so- 
lution is only known in the case of a piecewise linear 
potential [3]. In this article, we will derive the analytic 
tunneling solution for the piecewise linear and quadratic 
potential 
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(see Fig [I]) which does not make use of the thin- wall ap- 
proximation, see also recent work by [3]. However, one 
(non-differential) equation remains to be solved either 
numerically or in the limits of a large and small bubble 
radius. We show that our general analytical results re- 
duce to the thin-wall results and the results of [3] in the 
appropriate limits. For the most part of this work, we ex- 
clude the effects of gravity and thus ignore the magnitude 



of the true vacuum energy V—. We will only offer some 
qualitative arguments about the inclusion of gravity in 
the end, except for the thin- wall limit. 

We expect this tunneling solution to be the simplest 
yet realistic approximation to an arbitrarily shaped bar- 
rier exiting smoothly into a true minimum. The rea- 
son is that a quadratic potential towards the true vac- 
uum allows us to attach to its critical point smoothly 
a given shallowly sloped region of scalar potential con- 
taining a minimum. This may be particularly relevant 
for studying the dynamics of populating a landscape of 
slow-roll inflationary slopes via tunneling from some false 
vacuum. Any discussion of the relative prevalence of dif- 
ferent classes of inflationary models in a candidate funda- 
mental theory such as string theory will necessarily have 
to include the discussion of the bias incurred by popula- 
tion via tunneling. 

In the absence of gravity, the Euclidean action for a 
single scalar field with potential energy V{<p) is given by 
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and the bounce solution 4>(r) is determined by solving 
the Euclidean 0(4) symmetric equations of motion 
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<j>"{r) + -4>'{r) = d+V.. 



(3) 



where <f>' = d r (f> and the potential for our case is given by 
Eq. 0. 

Initially, the field is sitting in the false vacuum at 
0+ < 0, kept in place by a linear potential piece (shown 
as a dashed line in Fig [I]) attached to the left of the 
false minimum. Its slope is irrelevant in the following 
computation as long as it classically stabilizes the field 
in the false minimum. To start with, we first depict 
the numerically solved bounce solution schematically in 
Fig [2] The field sits at some value </>o (not necessarily 
close to the true vacuum 0_ = 1 in this case but to the 
right of the maximum of the potential) inside the bub- 
ble at r = 0. At the bubble radius Rt, the field crosses 
through the maximum of the potential at <P(Rt) = 0. 
Well outside the bubble, the field sits in the false vac- 
uum <j>(R+) = 0+ = -0.0001. 

The solutions 4>l and cf>R in the left and right part 
of the potential, respectively, have to fulfill the follow- 
ing boundary conditions: The bubble nucleates at rest 



2 




4> + 0_ 

FIG. 1: Schematic shape of the potential Eq. |TJ) 
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FIG. 2: The bounce solution for a = 10~ S ,A = 10~ J ,m = 
10 -6 , (f>- = l.The field sits at some value <f> = <f>o = 0.0037 in- 
side the bubble at r = 0. At the bubble radius, the field crosses 
through the maximum locus of the potential at (/>(Rt) = 0. 
Well outside the bubble, the field sits in the false vacuum 
<f>(R+) = <j>+ = —0.0001. For the definitions of a, A see the 
text below Eq. Q. 

<M0) = 0o > 0,<P'r(0) = 0. Outside of the bubble wall 
at finite r > R + , the field sits at rest in the false vacuum 
i.e 4> L (R+) = <j>+,<j>' L (R + ) = 0. Solving Eq. ^ with 
these boundary conditions for both parts of the poten- 
tial, we find 
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where I\(z) is the modified Bessel function of the first 
kind (see also equation (3.14) in pQ). 

In order to determine the constants <po, R+, and Rt, we 
use that the field configurations need to match smoothly 
at the bubble radius Rt with c/)l(Rt) = 4>b,(Rt) = and 
4>' l (Rt) = 4>' r (Rt)- Thus we need to solve 
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Calculating the Euclidean action Eq. ^ for the solutions 
Eq. Q, we can express both <j)Q and i?+ in terms of Rt 



to obtain 
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where we have introduced a — —<j) + /(j)^ > 0, and 
A = (—2\ + 4> + )/(m 2 4> 2 _) as a measure of the height of 
the potential barrier with values < A < 1. 
In order to find Rt, we combine Eq. ^ to get 
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which can be solved numerically. However, it is much 
more instructive to examine it in the limits of large and 
small ttiRt- But first we should briefly note that Eq. J8| 
can be used to remove the Bessel functions from Eq. (CT. 
We emphasize that the above expression for B now de- 
pends on the potential parameters and the only unknown 
quantity Rt, which is fully determined in terms of a, A 
by an algebraic Eq. This is one of our central results. 

Now we turn towards the two approximate solutions of 
our results. We start by taking the limit itlRt 3> 1, and 
in this limit, Eq. Q can be solved as 

R 3 + 40 1 ^ 

We note that for all allowed values of a and A, tuRt > § • 
In particular, we find that itlRt can be large either for 
the thin- wall limit A ps 1, or A < 1, but a > A (see 
Eq. (p0|). For the latter case, the potential minima are 
separated by sizable distances in the field space, as well 
as the potential energy between the false and true vac- 
uum. Therefore, the large tuRt limit encompasses more 
than the just thin-wall solution, see Table [I] For small 
differences in the vacuum energy e = (1 — /±)\m 2 (f) 2 _, in 
leading order of e, 
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which we will find to be identical to the result from using 
the thin-wall formalism. This is in accord with expecta- 
tions from the thin- wall formalism: The radius of the nu- 
cleating bubble grows toward oo as the minima become 
degenerate. Also, looking at the limit of the matching 
condition 0^ ; = in Eq. ^ 

4>o = <P- (l - sj\{mR T f /2 e- mR A , (11) 

it is clear that the true vacuum bubble nucleates close to 
the minimum of the quadratic part of the potential, <p ~ 
ef>—, for large values of tuRt- Plugging in the expression 
for Rt from Eq. ^ yields an unwieldy result for B. Thus 
we only display the expression for B in the e — > limit 
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Finally, we give an expression of B in terms of (f>o 
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In summary, in this limit, when either the thin- wall ap- 
proximations are valid or the minima are separated far 
away from each other, a large bubble nucleates close to 
the true vacuum. 

For the sake of comparison, we shall now give the es- 
sential results of the thin-wall calculation, where the nu- 



cleation rate is given by B 



2 7^ S f 
2e 3 



and for the potential 



Eq. ([!]) we can compute S\ directly, and in the small e 
limit it is 
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Plugging this into the expression for B we find the same 
expression as Eq. ( 12 ). The same way, we get Rt — 3Si/e 
which, after plugging in S\ and expanding in e, agrees 
with fit)]). 



We will now proceed to the opposite case of taking 
the limit of small mRx < 1. To second order in mRx, 
Eq. (|8| is solved by 
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where we discarded the negative solution. Now tuRt < 
[it can be equivalently casted as (0 < a < fi T /8) A (0 < 
A < {pip — 8a) 2 /16/x|, . Therefore small mRr limit cor- 
responds to the small values of a and A, where potential 
minima are closely spaced in the field space, but the po- 
tential difference between the false and the true minima 
is considerably large. In this limit, the exit point of the 
bubble can be conveniently expressed as 



t> = <t>- 1 
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and it shows that the bubble nucleates close to the tip 
of potential barrier which is intuitively understandable 
from the form the potential in this particular limit. We 
can express B for small tuRt in terms of a, A as 
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Again, we give an expression of B in terms of 4>q 
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In the limit where the quadratic part of the potential 
can be regarded as almost flat, the solution we found for 
small tuRt should agree with the results of a piecewise 



mRr B exact Bjjj Bthin-waii 
a = 0.01 0.6 0.0023 0.0022 72.4 
a = 0.1 1.2 0.3 0.04 113.3 
a = 0.5 2.5 23.9 0.6 529.8 

TABLE I: The mismatch of action: the linearized result Bdj 
of [3j compared to the thin-wall result Bthin-waii and the 
exact result found here. B is given in units of 4> 2 _ /m 2 and 
all values are quoted for A = 0.01. Clearly, there are regimes 
with tuRt = 0(1) where neither the linearized treatment 
of [3] nor the thin-wall approximation are sufficient. 



linear potential as studied by [3] . This is the case for the 
limit Ac 1. Taking the slope of the quadratic potential 
to be A_ = m 2 0_ , the parameter c defined by [3] becomes 
c = la I A. The value of the tunneling radius tuRt found 
by [3] is given by 
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in agreement with the corresponding expansion of 
Eq. (15). Expanding the expression for B in Eq. (|7|), 



we find perfect agreement with B from [3]. 

Let us recall what has been done so far. We reduced 
the computation of the tunneling amplitude for tunnel- 
ing from the linear to the quadratic part of the potential 
Eq. ([I]) to solving an exact Eq. pi . We provided approx- 
imate solutions for the limits of large and small tunnel 

L S 

radius mR T ' . Surprisingly, the smaller of the relative 
error between either the numerical result R^ and R^ or 
R T is always smaller than 10%, see Fig. [3] (a). This then 
defines globally the approximate solution tiiRt 



mR,T = < 



2VA 1 



A 



, A < (0.8a -0.5) 5 



(20) 



3+4a 1 
2 1-VA 



(0.8a -0.5) 2 < A 



In other words, we succeeded in computing the tunneling 
amplitude analytically in a non thin-wall solution with 
error better than 50%, see Fig. [3] (b). 

Finally we shall comment on the effects of gravity on 
the tunneling process. In the thin-wall limit [2] studied 
the case where either true or false vacuum are at zero 
energy, whereas [5] examined this issue for arbitrary val- 
ues of the false and true vacuum energy. Among other 
things, their work suggests the following intuitive idea. 
When the true vacuum is also in de-Sitter space, due to 
the "pull" on the bubble both from the inside and the 
outside of the wall, the nucleation rate should be en- 
hanced compared to the flat space limit. We can see this 
explicitly in the thin- wall approximation, i. e. large vtlRt 
and 1 — A C 1, where the effect of gravity is well under- 
stood. Following [5], we compute the effect of gravity on 
the exponent B of the tunnel rate, see Fig. [3] (c) . 

For sufficiently sub-Planckian values of 0_ we always 
have a constant asymptotic value for B / B , where B 
is without gravity. In the case of small = m 2 /V_ 



FIG. 3: The smaller of the relative error between the numerical value of Rt (panel a) and B (panel b) and its small or large 
mR T limit, Min(|i?^ - R^\jR%, \R% ~ JR%\/B%). The relative error is always smaller than 10% for R T and 50% for B. c) 
B/Bq of the tunneling rate for A = 0.95 in the presence of gravity as a function of 0_, with r\- = m 2 /V_. 



(i.e a very flat quadratic part) the ratio B/Bq is very 
small, giving a large gravitational enhancement of tun- 
neling. Contrary, for a steep potential with large r\- the 
tunneling probability is less enhanced. Increasing a sup- 
presses B further and thus enhances tunneling. This fits 
with the notion, that the gravitational correction is more 
important the thicker the barrier is in the field space, 
since for a given <p_ the barrier thickness increases with 
increasing a. Note, that for a ^> 1 and rj- <C 1 there is 
also Hawking-Moss tunneling possible [SHI]. 

In this work we have discussed quantum tunneling in 
field theory in a piecewise linear and quadratic scalar po- 
tential with a a false vacuum and a true vacuum. Such a 
potential is arguably the most simple yet realistic approx- 
imation to an arbitrarily shaped barrier exiting smoothly 
into a true minimum. The reason is that a quadratic po- 
tential allows us to attach to its critical point smoothly 
a given shallowly sloped region of scalar potential con- 
taining a minimum. Our result gives the tunneling rate 
in this situation exactly. Further, in the appropriate lim- 
its our result reduces to either the thin-wall result or 
the known result for piecewise linear potentials by [3]. 
However, there is a large region of barrier shape parame- 
ter space where neither of them is a good approximation 
to the full solution given here. The inclusion of further 
effects from approximating a given potential to higher 



than quadratic order, as well as a detailed incorporation 
of gravity in the non-thin-wall regime we leave as a topic 
for future work. Let us note in passing that in the context 
of meta-stable supersymmetry breaking vacua in gauge 
theories [5], choosing Nf = 3N c /2 flavors can imply a 
barrier shape with a -C A ~ 2/3. For a = 0.1 this 
gives itiRt — 9 which shows the linearized result of [3] 
with Bjjj ~ 2.5 x 10 2 overestimating the tunneling rate 
by an order of magnitude compared to our exact result 
B ~ 1.4 x 10 3 . Finally, we expect one relevant appli- 
cation of our results to be the dynamics of populating 
a landscape of slow-roll inflationary slopes via tunneling 
from some false vacuum. For this purpose, we have given 
the instanton action for the tunneling also as a function 
of cj>o, as the exit from the tunneling <p = <po, (f>o = forms 
the initial conditions for slow-roll inflation. 

Note added: While this paper was being finished, we 
became aware of |10j . which overlaps with our results. 
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